function r=rungekutta4(t0, t1, y0, n)
	h=(t1-t0)/n;
    % n: cantidad de pasos
    % y0: valor inicial
    % x0,x1: intervalo de evaluacion
	ts=t0:h:t1;
	y = y0;
	for i=1:n
		t = ts(i);
        it=i-1;
		k1 = h * f(t, y);
		k2 = h * f(t + h/2, y + (h/2)*k1);
		k3 = h * f(t + h/2, y + (h/2)*k2);
		k4 = h * f(t + h, y + k3);
        fprintf('\n%2.0f%10.6f%10.6f\n',it,t,y);
		y = y + (k1 + 2*k2 + 2*k3 + k4)/6;
	end
	r = y;
    fprintf('\n El punto aproximado y(x1) es = %8.6f\n',y)
    plot(ts,r,'o'),... 
    title('Runge Kutta 4'),... 
    xlabel('Tiempo'),ylabel('y = f(t)'),grid

